Discrete surface soiitons in two dimensions 
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We investigate fundamental localized modes in 2D lattices with an edge (surface). Interaction 
with the edge expands the stability area for ordinary soiitons, and induces a difference between 
perpendicular and parallel dipoles; on the contrary, lattice vortices cannot exist too close to the 
border. Furthermore, we show analytically and numerically that the edge stabilizes a novel wave 
species, which is entirely unstable in the uniform lattice, namely, a "horseshoe" soliton, consisting 
of 3 sites. Unstable horseshoes transform themselves into a pair of ordinary soiitons. 
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I. INTRODUCTION AND THE MODEL 

Soiitons on surfaces of fluids Q, solids @, and plas- 
mas Q are classical objects of experimental and theo- 
retical studies of nonlinear science. Recently, a new im- 
plementation of surface solitary waves was proposed Q 
and experimentally created [f| in nonlinear optics, in the 
form of discrete localized pulses supported at the edge 
of a semi-infinite array of nonlinear waveguides. Two- 
component surface soiitons were analyzed too jg, and it 
was predicted that soiitons may be supported at an edge 
of a discrete chain by a nonlinear impurity Ml- Parallel 
to that, surface gap soiitons were predicted [8( and cre- 
ated in an experiment 9] at an edge of a waveguiding 
array built into in a self-defocusing continuous medium. 
Very recently, the experimental creation of discrete sur- 
face soiitons sup por ted by the quadratic nonlinearity was 
reported as well [lg. In these cases, the soiitons are one- 
dimensional (ID). In Ref. a 2D medium with sat- 
urable nonlinearity was considered, with an embedded 
square lattice, that has a jump at an internal interface; 
there, stable asymmetric vortex soiitons crossing the in- 
terface were predicted, as a generalization of discrete vor- 
tices on 2D lattices |l2| and vortex soiitons supported by 
optically induced lattices in photorefractive media 
Soiitons supported by a nonlinear defect at the edge of a 
2D lattice were also considered Q. 

The search for surface soiitons in lattice settings is a 
natural problem, as, in any experimental setup, the lat- 
tice inevitably has an edge. In this paper, we report new 
results for surface soiitons in semi-infinite 2D lattices. We 
will first consider straightforward generalizations of local- 
ized modes studied previously on uniform lattices, viz., 
fundamental soiitons and two types of dipoles, oriented 
perpendicular or parallel to the surface. Then, we will 
introduce a novel species, the so-called horseshoe soliton, 
in the form of an arc abutting upon the lattice's edge. 
The existence, and especially the stability, of such a lo- 
calized mode is a nontrivial issue, as attempts to find 



a "horseshoe" in continuum media with imprinted lat- 
tices and an internal interface (similar to the medium in 
Ref. ^l|) have produced negative results [l5|. We find 
that, in the semi-infinite discrete medium, a horseshoe 
not only exists near the lattice edge, but perhaps more 
importantly has a stability region. For comparison, we 
also construct a family of the same wave pattern in the 
uniform lattice [which is a form of a stationary localized 
solution of the 2D discrete nonlinear Schrddinger (DNLS) 
equation that has not been discussed previously and we 
illustrate how/why it is of interest in its own right]. In 
particular, we find that this family of solutions is com- 
pletely unstable in the bulk of the uniform lattice, which 
stresses the nontrivial character of the surface-trapped 
horseshoes, in that they may be stabilized (for an appro- 
priate parameter range) by the lattice edge. 

The model of a semi-infinite 2D array of waveguides 
with a horizontal edge, that we consider below, is based 
on the DNLS equation for wave amplitudes u m n (z) in 
the guiding cores, z being the propagation distance: 

~t~ ^m.n— 1 4li m _ n ) -f- |tijn,n| ^m,n — 0; (1) 



for n > 2 and all to, where the prime stands for d/dz, 
and C is the coupling constant. At the surface row, 
n = 1, Eq. is modified by dropping the fourth 
term in the above parenthesis of Eq. JTJ (cf. the ID 
model in Refs. |f|), namely, w m .o = as there are no 
waveguides at n < 0. Note that, despite the pres- 
ence of the edge, Eq. admits the usual Hamiltonian 
representation, and conserves the total power (norm) , 

Stationary solutions to Eq. £Q| will be looked for as 
Mm, n = e lkz v m ^m where the wavenumber k may be scaled 
to 1, once C is an arbitrary parameter, and the stationary 
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solution obeys the equation 

+ Wm+l.n + V m -\,n - 4u m , n ) = 0, (2) 

with the same modification as above at n = 1. 

We will first report results of an analytical approxima- 
tion for the shape and stability of dipoles and "horse- 
shoes", valid for a weakly coupled lattice (C — > 0). This 
will be followed by presentation of corresponding numeri- 
cal results. Finally, we will briefly discuss the interaction 
of vortices with the lattice's edge. 

II. PERTURBATION ANALYSIS 

Analytical results can be obtained for small C, starting 
from the anti-continuum (AC) limit, C = (see Ref. fig 
and references therein). In this case, solutions to Eq. @ 
may be constructed as a perturbative expansion 

oc 

V m ,n — u v m,n- 
k=0 

In the AC limit proper, the seed solution, Vm] n , is zero 
except at a few excited sites, which determine the config- 
uration. 

By means of the analytical method, we will consider 
the following configurations: (a) a fundamental surface 

soliton, seeded by a single excited site, vf l = 1 (the first 
subscript 1 denotes the soliton's location in the horizon- 
tal direction), (b) surface dipoles, oriented perpendicular 
(bl) or parallel (b2) to the edge, each seeded at two sites, 

KU°0={-M},or KK1}={-1,1}, (3) 

and (c) the "horseshoe" 3-site structure 

{v^A%v%}={e m -,e^,e^}, (4) 

with #1,1 = 0, 0o,2 = 7T) 0-1,0 = 27T. We note in pass- 
ing that stable dipole states on the infinite lattice were 
predicted in Ref. [Ig], and later observed experimentally 
in a photorefractive crystal All the above seed con- 
figurations are real, and the horseshoe may, in principle, 
also be regarded as a truncated quadrupole, which is a 
real solution as well [l^ . 

At small C > 0, it is straightforward to calculate cor- 
rections to the stationary states at the zeroth and first 
order in C. Then, the stability of each state is deter- 
mined by a set of eigenvalues, A, which are expressed in 
terms of eigenvalues fi of the Jacobian matrix, to be de- 
rived in a perturbative form, M = J^'kLoC'Mk, from 
the linearized equations for small perturbations around 
a given stationary state |l6j . The stability condition is 
3?( A) = for all A, for the Hamiltonian system of interest 



herein (since if A is an eigenvalue, so are —A, A* and — A*, 
where the asterisk denotes complex conjugation). 

For the dipole and horseshoe configurations, (bl,b2) 
and (c), the calculations result in 

mo = c(~\ _\\ +o(c 2 ), 



M {c) = C 2 2-1-1 +0(C 3 ) (5) 

V 2 1 - i ) 

(the matrices for (bl) and (b2) coincide, at this order). 
From here, we obtain stable eigenvalues, A^ = 0, A 2 = 
±2y/Ci + 0(C), and = 0, \ { 2 c) = O (C 2 ) , A^ c) = 
±2V3Ci + (C 2 ). Both for the dipoles and half-vortex, 
one eigenvalue is exactly zero, as this corresponds to the 
Goldstone mode generated by the phase invariance of the 
underlying DNLS equation. As for eigenvalue A 2 > h be- 
comes different from zero at order O (C 2 ), and, as shown 
below, it plays a critical role in determining the stability 
of the horseshoe structure. 

We do not consider here the fundamental soliton, (a), 
as its destabilization mechanism is different (and requires 
a different analysis) from that of the dipoles and horse- 
shoes; in particular, the critical eigenvalues bifurcate not 
from zero, but from the edge of continuous spectrum (see 
below) . 

III. NUMERICAL RESULTS 

To examine the existence and stability of the above 
configurations numerically, we start with the fundamen- 
tal onsite soliton at the surface, (a). Basic results for 
this state are displayed in Fig. ^ At C = 0, there is a 
double zero eigenvalue due to the phase invariance. For 
small C > 0, this is the only eigenvalue of the lineariza- 
tion near the origin of the spectral plane (5ft(A), 3(A)). 
As C increases, one encounters a critical value, at which 
an additional (but still marginally stable) eigenvalue bi- 
furcates from the edge of the continuous spectrum, as 
mentioned above. With the further increase of C, this 
bifurcating eigenvalue arrives at the origin of the spectral 
plane, and subsequently gives rise to an unstable eigen- 
value pair, with Si(A) ^ 0, see Fig. ^ This happens for 
C > 1.41; we note in passing that the results reported 
herein have been obtained for lattices of size 10 x 10, but 
it has been verified that a similar phenomenology persists 
for larger lattices of up to 25 x 25. For comparison, we 
also display, by a dashed-dotted line, the critical unsta- 
ble eigenvalue for a fundamental soliton on the uniform 
lattice (as a matter of fact, for a soliton sitting far from 
the edge), which demonstrates that the interaction with 
the edge leads to conspicuous expansion of the stability 
interval of the fundamental soliton. This may also be 
justified intuitively, as the instability of the fundamen- 
tal soliton emerges closer to the continuum limit which 
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FIG. 1: (Color online) Dynamical features of the fundamental 
surface soliton. The top left and right panels show, respec- 
tively, the soliton's norm, P, and the real part of the critical 
stability eigenvalue, versus the lattice coupling, C. For com- 
parison, the dash-dotted lines show respective quantities for 
a fundamental soliton in the uniform lattice. The instability 
is due to an eigenvalue pair bifurcating from the edge of the 
phonon band, that eventually crosses the origin of the spec- 
tral plane and becomes real. In the bottom panels, we display 
the linear-stability spectrum of a single-site soliton, and snap- 
shots of its evolution (contour plots of \u m ,n\ 2 ), for C = 1.43, 
slightly above the instability threshold. 

is well-known to be unstable to (very slow) collapse type 
phenomena. A surface variant of the relevant structure 
enjoys the company of fewer neighboring sites and hence 
is "more slowly" (in parameter space) approaching its 
continuum limit, in comparison with its infinite lattice 
sibling. 

Development of the instability of the fundamental sur- 
face soliton (in the case when it is unstable) was exam- 
ined in direct simulations of Eq. Q). As seen in Fig.^ hi 
this case the excitation propagates away from the edge, 
expanding into an apparently disordered state. This is 
justifiable as for this parametric regime there is no stable 
localized state neither in the vicinity of the surface or in 
the (more unstable) bulk of the uniform lattice. 

Next, in Fig. [21 we present results for the vertical and 
horizontal dipoles, (bl) and (b2) (see top left and mid- 
dle panels of Fig. [21 seeded as per Eq. J3J). At C = 0, 
the dipole has two pairs of zero eigenvalues, one of which 
becomes finite (remaining stable) for C > 0, as shown 
above in the analytical form. Our numerical findings re- 
veal that, in compliance with the analytical results, the 
dipoles of both types give rise to virtually identical finite 
eigenvalues (hence only one eigenvalue line is seen in the 
left middle panel of Fig- - As seen in the right middle 
panel in Fig. [21 both dipoles lose their stability simulta- 
neously, at C « 0.15. Continuing the computations past 
this point, we conclude that the vertical and horizon- 
tal dipoles become different when C attains values ~ 1. 




FIG. 2: (Color online) Top panels: Vertical (bl) and horizon- 
tal (b2) dipoles, which, at C = 0, are seeded through Eq. JHJ. 
The top left and middle panels show examples of these dipoles 
at C = 1, and the right panel shows their norms. Middle and 
bottom panels: The same as Fig.Q but for the configurations 
bl and b2. The vertical dipole, bl, disappears via a saddle- 
node bifurcation at C ~ 2.17. The left middle panel displays 
the eigenvalue bifurcating from zero at C = 0, the dashed 
line being the analytical approximation described in the text, 
i.e., S5(A) = 2yC. The right middle panel shows the onset of 
the bl (solid lines) and b2 (dashed lines) instabilities for the 
dipole modes, as found from numerical computations. An ex- 
ample of the spectrum and nonlinear evolution of an unstable 
vertical dipole are presented in the bottom panels, for C = 1. 

Eventually, the (already unstable) vertical configuration, 
bl, disappears in a saddle-node bifurcation at C w 2.17, 
while its horizontal counterpart, b2, persists through this 
point. Furthermore, there is a critical value of C at which 
an eigenvalue bifurcates from the edge of the continuous 
spectrum. Eventually, this bifurcating eigenvalue crosses 
the origin of the spectral plane, giving rise to an unstable 
eigenvalue pair, with 5R(A) ^ 0. The value of C at which 
this secondary instability sets in is essentially smaller for 
bl, i.e., C w 1.55, than C w 2.61 for b2. We thus con- 
clude that the horizontal dipole, b2, is, generally, more 
robust than its vertical counterpart, bl. This conclusion 
seems natural, as the proximity to the edge stabilizes 
the fundamental soliton (as shown above), and in the 
horizontal configuration the two sites that constitute the 
dipole are located closer to the border. 

Nonlinear evolution of unstable dipoles was examined 
too, in direct simulations. As seen in the example (for 
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FIG. 3: (Color online) The same as Fig.[5]for the "horseshoe" 
configuration, Eq. J3) . The top panels show an example of the 
structure and its norm as a function of C (black solid curve). 
The solid curves in the middle panels display the real and 
imaginary parts of the key stability eigenvalues (the dashed 
fine in the right middle panel is the analytical approximation 
for the imaginary part presented in the text). For comparison, 
the red dash-dotted lines show the same characteristics for a 
family of horseshoe solitons created in the uniform lattice. It 
is seen that the latter family is completely unstable, while the 
same waveform trapped at the edge of the lattice has a well- 
defined stability region. This is a remarkable consequence 
of the presence of the surface. The bottom panels present 
the linear instability spectrum of the horseshoe state with 
C — 0.26, and its evolution in time. 

configuration bl) shown in the bottom panels of Fig. [21 
the instability again results in a disordered state. 

Proceeding to the completely novel configuration (c) 
of the horseshoe, we note that, because it was seeded at 
three sites when C = 0, there are three pairs of zero eigen- 
values in the AC limit. Above, it was shown analytically 
that one pair of these eigenvalues becomes finite at order 
0(C), and another at 0(C 2 ). Numerical results demon- 
strate that the first pair remains stable until it collides 
with the edge of the continuous spectrum, which happens 
at C » 0.25. As mentioned above, the second eigenvalue 
pair, bifurcating from zero at order O (C 2 ) , is critical for 
the stability of configuration (c). The numerical results 
show that this pair bifurcates into a stable (i.e., imagi- 
nary) one, and, as shown in Fig. EI the horseshoe remains 
stable up to C ~ 0.25. 

To understand the effect of the surface on the stability 
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FIG. 4: (Color online) The supersymmetric vortex cell seeded 
as per Eq. Q. The top left and right panels show, respec- 
tively, the real and imaginary parts of the solution and its 
(in)stability spectrum for C = 0.4. The bottom panel dis- 
plays imaginary and real parts of the stability eigenvalues 
versus C. The solid and dashed lines show numerical and 
analytical results for small C. For comparison, red dashed- 
dotted lines depict the same numerically found characteristics 
for a supersymmetric vortex on the infinite lattice. 



of the horseshoes, it is relevant to compare them to their 
counterparts in the infinite lattice, i.e., similar structures 
that may be created, starting from the AC limit taken 
as per Eq. (@J, far from the domain edge. By itself, the 
latter is a novel family of localized solutions to the DNLS 
equation in 2D. However, the important feature is that, 
the structure in the bulk of the uniform lattice (unlike the 
quadrupole that may be stable), is unstable for the entire 
family of horseshoes [in fact, the 0(C 2 ) eigenvalue pair, 
bifurcating from zero at C = 0, immediately becomes real 
in this case, see the corresponding, quadratically growing 
dashed-dotted line in the left middle panel in Fig. 
Thus, the horseshoe attached to the surface is another 
example of the localized mode stabilized by the lattice's 
edge (however, unlike the above example of the stabilized 
fundamental soliton, the dependence of the horseshoe's 
stability on the border is crucial, as it may never be stable 
in the infinite lattice). 

The bottom panels of Fig. [3] exemplify the evolution of 
the horseshoe when it is unstable. Unlike the situations 
with configurations (a) and (b), the unstable horseshoe 
does not decay into a disordered state, but rather splits 
into a pair of two fundamental solitons (with S = 0) , one 
trapped at the surface and one found deeper inside the 
lattice. 
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IV. SURFACE EFFECT ON THE EXISTENCE 
OF VORTICES 

Besides the stabilization effect reported above, the lat- 
tice edge may exert a different effect, impeding the exis- 
tence of some solutions. As an example, we consider the, 
so-called supersymmetric [lfij , lattice vortex attached to 
the edge, i.e., one with the vorticity (S — 1) equal to 
the size of the square which seeds the vortex at C = 0, 
through the following set of four excited sites, cf. Eq. 

{v£l v ( °lA% «$} = {«^> ) e*.S^ I ^" } , (6) 

with 6*0,1 = 0, M = tt/2, 0i,2 = 7T, and O , 2 = 3tt/2 
(unlike the above configurations, this one is not a purely 
real one). While supersymmetric vortices exist in uni- 
form lattices (including anisotro pic ones) and have their 
own stability regions there [HI, l20j|. numerical analysis 
shows that the localized mode initiated as per Eq. in 
the model with the edge cannot be continued to C > 0. 
In fact, we have found that, to create such a state at fi- 
nite C, we need to seed it, at least, two sites away from 
the edge, i.e., as 

KWMM2} = {^/ l -',e^ ) e^} . (7) 

[the set on the right-hand side is a translated version of 
that of Eq. ©]. Numerically found stability eigenval- 
ues for this structure are presented in Fig. 01 along with 
the analytical approximation, obtained by means of the 
same method as above. In all, there are 4 pairs of analyti- 
cally predicted eigenvalues near the spectral plane origin 
(given the four initial seed sites of the configuration). 



More specifically, these are: A = (associated with the 
phase invariance), A = ±2Ci (a double eigenvalue pair), 
and A = ±\/32C 3 i (a higher order eigenvalue pair). As 
it can be seen in the figure, already at this distance of 
two sites from the boundary, the behavior is sufficiently 
close to that of an infinite lattice. 



V. CONCLUSION 

We have demonstrated that properties of fundamen- 
tal localized modes in the 2D lattice with an edge may 
be drastically different from well-known features in the 
uniform lattice. In particular, the edge helps to increase 
the stability region for the ordinary solitons, and induces 
a difference between dipoles oriented perpendicular and 
parallel to the lattice's border. Additionally, regular su- 
persymmetric vortices cannot be created too close to the 
border. Most essentially, the edge stabilizes a new species 
of solitons which is entirely unstable in the uniform lat- 
tice, the so-called "horseshoe solitons". In that sense, 
the presence of the surface may produce a very desirable 
"stability broadening" impact on a number of different 
solutions; this may become especially interesting in appli- 
cations related to waveguide arrays, among other optical 
or soft-matter systems. 

Natural issues for further consideration are horseshoes 
of a larger size, and counterparts of such localized modes 
in 3D lattices near the edge. In the 3D lattice, one can 
also consider solitons in the form of vortex rings or cubes 
[2l| set parallel to the border. 
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